A transition from river networks to scale-free networks 
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A spatial network is constructed on a two dimensional space where the nodes are geometrical 
points located at randomly distributed positions which are labeled sequentially in increasing order 
of one of their co-ordinates. Starting with N such points the network is grown by including them 
one by one according to the serial number into the growing network. The t-th point is attached to 
the i-th node of the network using the probability: 7Ti(t) ~ ki(t)l% where fei(t) is the degree of the 
i-th node and £u is the Euclidean distance between the points t and i. Here a is a continuously 
tunable parameter and while for a = one gets the simple Barabasi-Albert network, the case for 
a —* — oo corresponds to the spatially continuous version of the well known Scheidegger's river 
network problem. The modulating parameter a is tuned to study the transition between the two 
different critical behaviors at a specific value a c which we numerically estimate to be -2. 
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Scale-free networks (SFN) are highly inhomogeneous 
with a power law decay of their nodal degree distribu- 
tions signifying the absence of a characteristic value for 
the nodal degrees [lj]. Extensive research over last sev- 
eral years revealed that such networks indeed occur in 
different real-world systems like protein interaction net- 
works in Biology, Internet and World-wide web (WWW) 
in electronic communication systems, airport networks 
in public transport systems etc. @, y, 0| • On the other 
hand, river networks are relatively simple spatial net- 
works which were being studied over last several decades 
from the geological point of view. During the last decade 
or so physicists have also studied properties of river net- 
works with many different simple model networks mainly 
from the interests generated about their fractal proper- 
ties, a popular topic of critical phenomena 

In this paper we report our study of an weighted spatial 
network on the two-dimensional plane. The weight of a 
link is evidently the Euclidean length of the link. Tuning 
a parameter which modulates the strength of the contri- 
bution of the link weight in the attachment probability, 
we are able to obtain networks similar to the directed 
river network model in one limit. On the other limit of 
the parameter we obtain scale-free networks. The tran- 
sition point for the crossover between the two types of 
behaviors is studied. 

The simplest river network model on a lattice is the 
Scheidegger's river network with a directional bias 
This is simply described on an oriented square lattice: 
Each lattice site is associated with an outgoing arrow 
representing the direction of the flow vector from that 
site. Only two possible choices for this arrow are pos- 
sible: it may direct to the lower left lattice site or to 
the lower right. An independent and uncorrelated as- 
signment of an arrow from each site results a 'Directed 
Spanning Tree' (DST) network @, [!]• Such networks are 
characterized mainly by the critical exponents associated 
with the distributions of the river basin area as well as 
the length of the longest river at each site. The set of as- 



sociated exponents constitute the universality class of the 
Scheidegger's river network which are different from the 
similar exponents of the isotropic spanning tree networks 
1. 

On the other hand while studying the scale-free prop- 
erties of different real-world networks Barabasi and Al- 
bert (BA) argued that there is a 'Rich get richer' mech- 
anism in-built with the growth process of every SFN [1]. 
They proposed a model of generating scale-free networks 
where new nodes are introduced to the growing network 
at a rate of one per unit time step which are connected 
with the growing network with m distinct links with a 
probability proportional to the individual nodal degree: 
7Tj(i) oc fcj(i) [1]. Also there are some other directed scale- 
free networks whose links are meaningful only when there 
is a connection from one end to the other but not along 
the reverse direction, e.g., the World-wide web lldl . the 
phone-call network [llj and the citation network [12j etc. 

Real-world networks whose nodes are geographically 
located in different positions on a two-dimensional Eu- 
clidean space are very important in their own right. For 
example the electrical networks in power transmissions, 
railway or postal networks in transport networks are few 
of the very well known spatial networks. Research over 
last few years have also revealed that two very impor- 
tant spatial networks like the Internet [3, [3, EH and 
the Airport networks 1|| 17 1 have scale-free structures. 

Weighted networks are those whose links are associ- 
ated with non-uniform weights Wij . Therefore the spa- 
tial networks are by definition weighted networks whose 
link weights are the Euclidean lengths of the links. For a 
weighted network one can define the strength of a node 



as the total sum of the weights Sj 



YijWi 



How the 



average nodal strength (s(k)) varies with the degree k, 
i.e., (s(k)) oc k$ is also a non-trivial question. Non- 
linear strength-degree relations have been observed for 
the Internet as well as the Airport networks. In this con- 
text a detailed knowledge of link length distribution is 
also important e.g., in the study of Internet's topologi- 
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FIG. 1: Pictures of the network generated for the four dif- 
ferent values of the modulating parameter a using the same 
set of N = 1024 randomly distributed points within the unit 
square box. Link lengths are large for a — which gradually 
shrink in length as a decreases and the network becomes the 
directed minimal spanning tree at a — — oo. 
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FIG. 2: (Color online) The average length of a link (£(a, N)) 
has been plotted in (a) with the continuously tunable param- 
eter q for different values of the system sizes: N — 2 8 , 2 10 , 2 12 
and 2 14 . N increases from top to bottom. In (b) a scaling 
is shown with the system size dependent critical value of a c 
which approaches -2 as N — » oo. 



cal structure for designing efficient routing protocols and 
modeling Internet traffic. Early studies like the Waxman 
model describes the Internet with exponentially decaying 
link length distribution [18}. Yook et. al. observed that 
nodes of the router level network maps of North America 
are distributed on a fractal set and the link length distri- 
bution is inversely proportional to the link lengths [l5j |. 
A number of model networks on the Euclidean space have 
been studied in different contexts IjJ 20j . 

We consider here a stack of nodes dropped one by one 
on a substrate with increasing vertical co-ordinates. Each 
node is connected randomly with a specific link length 
dependent probability of attachment to a node of the 
already grown stack. 

A network of N nodes is grown within an unit square 
box on the two-dimensional x — y plane. Nodes are 
represented by N points selected at random positions 
(xi,yi),i — l,N by generating their values from an uni- 
form probability distribution {0, 1}. The first point is 
placed by hand at the bottom of the box with yi = 0. 
All other points are assigned serial numbers in increas- 
ing order of their y-coordinates: y\ < yi < 2/3... < yjy. 
We use the geometry of a cylinder i.e., impose the open 
boundary condition along the y-direction but the periodic 
boundary condition along the x-direction. This implies 
that the space is continuous along the x-direction and 
any node very close to the x = 1 line may have a right 
neighbor inside the box and very close to the line x = 
and vice-versa. 

To start with we assume that the first node at the 



bottom of the box has a 'pseudo' degree fci = 1. Then 
the nodes from 2 to N are connected to the network 
by one link each. The time t measures the growth of 
the network by the number of nodes. The t-th node is 
then linked to the growing network with an attachment 
probability 



7T<(t) 0C ki{t)ll 



(1) 



where In denotes the Euclidean distance between the t- 
th and the i-th nodes maintaining the periodic boundary 
condition. This implies that the attachment probability 
has two competing factors. The linear dependence on 
the degree ki(t) enhances the probability of connection 
to a higher degree node where as the factor £ a reduces 
the probability of selection when a > and enhances 
when a < 0. The special case of a = the attach- 
ment probability in Eqn. (1) clearly corresponds to the 
Barabasi-Albert model. We now discuss the properties 
of the network by continuously tuning the parameter a 
through its accessible range. 

In the limit of a — > —00 every node connects its nearest 
node i in the downward direction corresponding to the 
smallest value of the link length I with probability one 
irrespective of its degree ki and therefore the probability 
of attachment to any other node is identically zero. This 
link may be directed either to the left or to the right 
depending on the position of the nearest node (Fig. 1). 

Let us first study the first neighbor distance distri- 
bution. Consider an arbitrary point P at an arbitrary 
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FIG. 3: The degree probability distribution P(k) vs. k of 
the network in the limit of a — » — oo and for TV = 2 18 . The 
solid line is a fit and its functional form is given in Eqn. (7). 



position. The probability that its first neighbor is posi- 
tioned on the semi-annular ring within r and r + dr in 
the downward direction (which can be done in N — 1 dif- 
ferent ways) and all other N — 2 points are at distances 
larger than r + dr is: 



Prob(r, N)dr = (N - 1) urdr (1 - 7rr 2 /2) 



N-2 



(2) 



In the limit of N — > oo it can be approximated that 
/Y-1«JV-2«W and since the average area per 
point decreases as 1/N, irr 2 /2 is very small compared 
to 1. Therefore (1 — irr 2 /2) N ~ 2 is approximated as 
cxp(— irNr 2 /2). In the limit of N — ► oo the probabil- 
ity density distribution is therefore: 

Prob(r, N) = irNr exp(-7r7Vr 2 /2) (3) 

or in the scaling form: 

1 



Prob(r) = [irrVN] exp(-7r[rV/V]72) = £(x)(4) 



where the scaling function Q{x) = ttx exp(— irx 2 /2) and 
x = ryN. Numerical results for the link length distri- 
bution of different system sizes verifies this distribution 
very accurately. 

The typical length (£(a, N)) of a link for a network of 
size N and generated with a specific value of the param- 
eter a is estimated by averaging the link length over all 
N — 1 links of a network as well as over many indepen- 
dent configurations. For this purpose one can define a 
total cost function C(a,N) of the network which is the 
total length of all the links: 



C(a,N) 



Zfj^ij and (t(a, N)) 



C(a,N)/(N-l).(5) 



In Fig. 2(a) we show the variation of the average link 
length with a. The (£(a, N)) — > as a — > — oo and 
gradually increases with increasing a. Around a c = —2 
the (£(a, N)) increases very fast and finally approaches 
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FIG. 4: (a) Variation of the scaling exponents Vk(a), /ifc(a) 
and 7fc(a) characterizing the degree distribution in Eqn. (8) 
with the modulating parameter a. (b) Variation of the expo- 
nents T a (a), Ti(a) and /3(a) with a. 



unity as a — ► +oo. The steep growth around a c be- 
comes increasingly sharper with increasing system size. 
For very large networks N — > oo it appears that for 
a < a c (£(a, N)) — > and for a > a c it approaches 
a finite value. Such a system size dependence is quanti- 
fied by a finite-size scaling of this plot as shown in Fig. 
2(b). The data collapse shows that 



{£{a, AO) oc N-° m G[{a - a c (N))N° 



(6) 



The critical values of a c (N) for a system of size iV is lo- 
cated at the value of a where (£(a, TV)) increased most 
rapidly with a. The values of a c (N) so obtained are 
-1.37, -1.46. -1.54 and -1.60 for TV = 2 8 , 2 10 , 2 12 and 2 14 
and are extrapolated with iV _al to get a c = a c (oo) = 
—2.0 ± 0.10. Similar results of a c = —2 have been ob- 
tained in [UIH]. 

In the limit of a — > — oo and N — > oo, the fractions 
of nodes with degrees 1, 2 and 3 are found to be 0.213, 
0.586 and 0.192 respectively and decreases very fast for 
higher degree values. The whole distribution fits nicely 
to a sharp Gamma distribution as (Fig. 3): 



P(k) oc k' b exp(-4:.2k). 



(7) 



In the range a > a c we observed that the network has 
a scale-free structure. For large system sizes the degree 
distribution follows a power law like P(k) oc fc -7 but for 
finite systems a finite-size scaling seems to work well: 



Prob(/c,iV) oc N-^Tuik/N 1 - 



(8) 
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FIG. 5: (Color online) Plot of the river basin area a distri- 
bution in our model for different network sizes N . In (a) 
Prob(a, N) has been shown for three different values of a: 
a = -co at the top, a — a c at the middle and a — at the 
bottom. In (b) a finite size scaling of Prob(a, N) with N has 
been shown for a = a c for N = 2 s , 2 9 and 2 10 . 



The scaling function Fk{x) ~ x~ lk for x — > and de- 
creases faster than a power law for a; >> 1 so that, 
7fe = Mfe/^fe- For a range of a values the scaling ex- 
ponents are measured and it is observed that all three 
exponents 7fe , fik and v% are dependent on the value of a 
(Fig. 4(a)). 

For the river network problems the size of the drainage 
area is a popular quantity to measure. The amount of 
water that flows out of a node of the river network is pro- 
portional to the area whose water is drained out through 
this node. On a tree network the drainage area etj is 
defined at every node i and is measured by the num- 
ber of nodes supported by i on the tree network. A well 
known recursion relation for a* is: a; = TljOjSij + 1 where 
the dummy index j runs over the neighboring nodes of 
i and Sij = 1 if the flow direction is from j to i, oth- 
erwise it is zero. The probability distribution Prob(a) 
of the drainage areas is the probability that a randomly 
selected node has the area value a. ft is known that for 
river networks this distribution has a power law variation: 
Prob(a) - a" r " [|. 

The drainage area distribution is measured first in the 
limit of a — ► — oo for our networks of different sizes N — 
2 12 ,2 14 and 2 16 . Direct measurement of the slopes of 
double logarithmic plots of Prob(a, N) vs. a gives values 
of the exponent r rs 1.33 which varied little with the 
system size (Fig. 5(a)). This estimate is consistent with 
that obtained from the finite size scaling analysis. An 
excellent scaling of the data over different system sizes is 



obtained as: 

Prob(a, N) cx N~^ a T(a/N Ua ) (9) 

and the exponent t q = ji a /v a . In the limit of a — > — oo we 
estimated fi a (— oo) » 1.00 and v a {— oo) w 0.75 giving a 
value for the exponent T a (— oo) rs 1.33. These values are 
very much consistent with the same exponents of Schei- 
degger's river network model where r a = 4/3 is known 
exactly Q. Similarly for a = we could reproduce the 
known values of // a (0) ~ 2.00 and v a {Q) « 1.00 with 
t<x(0) ~ 2.00 as obtained in [12} . Finally we measured 
the same distribution at the transition point a c = — 2 
and obtained /i Q (— 2) » 0.45 and v a {2) w 0.29 with 
r Q (-2) « 1.50 (Fig. 5(b)). 

Another quantity of interest is the length Li of the 
longest up-stream meeting at the node i. It's magni- 
tude is the number of links on the longest path ter- 
minating at i. Prob(L) therefore denotes the proba- 
bility that an arbitrarily selected node i has Lj = L. 
Given a network of size AT, L values are measured at 
every node and then the data is sampled over many 
uncorrelated network configurations and for different 
network sizes. A similar finite size scaling form like 
Prob(i,A) oc N-» L Q{L/N VL ) works here as well. We 
obtain ^(—aa) w 0.75, vl(— oo) w 0.50, tx(— oo) = 1.5 
and n L (2) w 0.2, v L {2) » 0.1, r L (2) a 2.0. 

Finally, we studied the strength-degree relation in our 
model. Here the weight of a link is the length of the link 
lij and therefore the strength of the node i is s,; = Sj^y. 
The strength (s(k)) per node averaged over all nodes of 
degree k of the network as well as over many independent 
realizations varies with the degree k as: (s(k,a)) oc k^^ . 
Numerically we observe that the exponent (3(a) varies 
with the tuning parameter a. Fig. 4(b) summarizes the 
variation of the three exponents T a (a), T£(a) and (3(a) 
within the range of a varying between -10 and 0. For a < 
—5, T a (a) and tx(oj) values coincide with their values at 
a = — oo. Between —5 < a < —2, tl(ol) slowly increases 
to 2 and T a (a) increases to 1.5. 

To conclude, we have defined and studied a network 
embedded in the Euclidean space. Random distribu- 
tion of nodes are sequentially numbered in increasing 
heights and the degree dependent attachment probability 
is modulated by the a-th power of the link length. This 
continuously tunable parameter a interpolates between 
the Scheidegger's river network and the Barabasi- Albert 
Scale-free network. It appears that there exists a critical 
value a c such that for a < a c the critical behavior of the 
network is like the Scheidegger's river network, whereas 
for a > a c critical exponents are indistinguishable from 
those of ordinary BA network. Our numerical study in- 
dicates a c is likely to be -2. 

Discussion with G. Mukherjee and K. Bhattacharya 
are thankfully acknowledged. 
Electronic Address: manna@bose.res.in 



[11 

[1] A.-L. Barabasi and R. Albert, Science, 286, 509 (1999). [12 

[2] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 

(2002). [13 

[3] S. N. Dorogovtsev and J. F. F. Mendes, Evolution of Net- 
works, Oxford University Press, 2003; M. E. J. Newman, [14 
SIAM Review 45, 167 (2003). 

[4] R. Pastor-Satorras and A. Vespignani, Evolution and [15' 
Structure of the Internet, A Statistical Physics Approach, 
Cambridge University Press, 2004. [ig 

[5] I. Rodrguez-Iturbe, A. Rinaldo, Fractal River Basins: 

Chance and Self- Organization, Cambridge University [17 
Press, 2001. 

[6] A.E. Scheidegger, Geol. Soc. Am. Bull. 72, 37 (1961); [ig 

Water Resour. Res. 4, 167 (1968). 
[7] F. Harary, Graph Theory, Addison- Wesley Publishing [19 

Company, Inc., Reading, Mass., 1969. 
[8] A. G. Bhatt and R. Roy, Appl. Probab. 36, 19 (2004). [20 
[9] S. S. Manna, D. Dhar and S. N. Majumdar, Phys. Rev. [21 

A. 46 R4471 (1992). [22 
[10] S. Lawrence and C. L. Giles, Science, 280, 98 (1998); 

Nature, 400, 107 (1999), R. Albert, H. Jeong and A.-L. 



Barabasi, Nature, 401, 130 (1999). 

W. Aiello, F. Chung and L. Lu in Proc. 32-nd ACM 
Symp. Theor. Comp. (2000). 

P. L. Krapivsky and S. Redner, Phys. Rev. E. 63, 066123 
(2001). 

M. Faloutsos, P. Faloutsos and C. Faloutsos, Proc. ACM 
SIGCOMM, Comput. Commun. Rev., 29, 251 (1999). 
R. Pastor-Satorras, A. Vazquez and A. Vespignani, Phys. 
Rev. Lett. 87, 258701 (2001). 

S. H. Yook, H. Jeong and A.-L. Barabasi, Proc. Natl. 
Acad. Sci. (USA) 99, 13382 (2002). 
R. Guimera and L. A. N. Amaral, Eur. Phys. Jour. B, 
38, 381 (2004). 

A. Barrat, M. Barthelemy, R. Pastor-Satorras, A. Vespig- 
nani, Proc. Natl. Acad. Sci. (USA), 101, 3747 (2004). 

B. Waxman, IEEE J. Selec. Areas Commun., SAC, 6, 
1617 (1988). 

A. F. Rozenfeld, R. Cohen, D. b-Avraham and S. Havlin, 
Phys. Rev. Lett. 89, 218701 (2002). 

S. S. Manna and P. Sen, Phys. Rev. E 66, 066114 (2002). 
J. M. Kleinberg, Nature, 406, 845 (2000). 
M. R. Roberson and D. ben-Avraham, Phys. Rev. E 74, 
017101 (2006). 



